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In the example of the Schrodinger/KdV equation we treat the theory as equivalence 
of two concepts of Liouvillian integrability: quadrature integrability of linear differ- 
ential equations with a parameter (spectral problem) and Liouville's integrability 
of finite-dimensional Hamiltonian systems (stationary KdV-equations) . Three key 
objects in this field: new explicit ^-function, trace formula and the Jacobi problem 
provide a complete solution. The 0-function language is derivable from these ob- 
jects and used for ultimate representation of a solution to the inversion problem. 
Relations with non-integrable equations are discussed also. 
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1. Introduction 

As it is understood in contemporary language, the theory of algebraic-geometric 
(or finite-gap) integration is a theory of integration of integrable nonlinear (1+1)- 
soliton partial differential equations (pde's), finite-dimensional Hamiltonian dy- 
namical systems, and spectral problems defined by ordinary differential equations 
(ode's). A large number of papers devoted to these problems, appeared over the 
last three decades, show that the vague term 'integration of integrable' does not 
mean an automatical integrating procedure in some simple sense of the word. The 
contributions that are still being made to this theory testify to its vitality: all the 
evidence points to the continuance of its growth. 

Although the term 'completely integrable' is quite appropriate one, there is not 
commonly agreed answer to the questions about integrability and, in particular, 
to the question in the title of the present work. For example, in the theory of 
dynamical systems this means the well defined Liouvillian integrability, but the 
search for separability variables (followed by action-angle ones) is a subject of an 
independent theory (Flaschka & McLaughlin 1976; Vanhaecke 1996). In the case 
of (l+l)-integrable pde's, the theory is, in fact, a treatment of these equations as 
infinite-dimensional analogues of Hamiltonian systems (Zakharov & Faddeev 1971; 
Gardner 1971) plus set of nontrivial exact solutions in terms of elementary functions: 
solitons, positons, and their relatives (Ablowitz & Segur 1981). There are natural 
multidimensional generalizations of the theory (the Kadomtsev-Petviashvili (kp), 
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Davey-Stewartson equations, their hierarchies, etc.) which are rather well devel- 
oped also. Integrability of spectral problems (direct and inverse) is usually associ- 
ated with the inverse scattering method, soliton theory (Ablowitz & Segur 1981; 
Novikov 1984; Levitan 1987) and its 9-generalizations (Dubrovin 1975; Matveev 
1976; McKean & van Moerbeke 1975; Flaschka 1975; Novikov 1984; Belokolos et all 
1994; Gesztesy & Holden 2003). We must include here the deep links of this field 
to diverse areas of mathematics shedding light on the mechanism of integrabil- 
ity: Hirota's r- function method, Painleve analysis (Conte 1999), Abclian varieties 
(Dubrovin 1976; Vanhaecke 1996), Darboux-Backlund algebraic transformations 
(Matveev & Salle 1991), differential geometry (Darboux 1915), theory of commut- 
ing differential operators (Krichever 1978) etc. Each of these links is sufficient in 
itself to provide a complete development; combined, they exhibit an unusual wealth 
of ideas and furnish rich resources of new interrelations. 

Nevertheless, there is a common object for all the approaches mentioned above: 
a 0-representation for a wide class of solutions and fundamental ^-function as 
a solution to associated linear spectral problem. The key role of this object was 
discovered by Its & Matveev (1975) in the example of the KdV-equation. The 
next years after, Krichever (1976, 1977a) put such a construction into a basis of 
integration of soliton equations and it became clear that this link between KdV- 
equation and the Schrodinger operator 

*"-u* = A*, u = <j){x) (1.1) 

is not an exception but a common feature. The ideology is spread to the whole 
class of such problems. Presently this is known as a concept of the Baker-Akhiczer 
function (Krichever 1978) and such an approach generates all the (l+l)-soliton 
equations and their hierarchies if an algebraic curve has been specified. 

Originally the term 'finite-gap' meant spectral problem for the smooth real 
periodic potential with finitely many lacunae at spectrum of Schrodinger's operator. 
Later, in the 1970's, in works by Matveev, Its, Dubrovin, Krichever, and others that 
treatment was generalized to quasi-periodic complex valued potentials and related 
to methods of algebraic geometry and 0-functions (Matveev 1976; Dubrovin et all 
1976; Krichever 1977b). By this reason throughout the paper we keep traditional 
terminology 'finite-gap' but identify it with 'algebraic-geometric'. 

It is difficult to keep pace with the continuing growth of the literature which is 
due to the activity of mathematicians. To become acquainted with the background 
of this field, as well as find a complete set of references, it is perhaps best to consult 
the surveys written by initiators of the theory. 

2. Solvable potentials for the Schrodinger equation 

Until the pioneer work by Novikov (1974) the Schrodinger equation (1.1) was mostly 
an object of the spectral theory of operators (Akhiezer 1961; Levitan & Sargsjan 
1975). By that time few solvable (in different senses) examples were known. 

1. The constant potential (trivial case); 

2. Quantum harmonic oscillator u = x 2 ; 

3. Linear potential u — x; 
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4. Decaying potentials u = —n (n+ 1) cosh 2 x and their rational degenerations; 

5. Reflectionless potentials of Bargmann (solitons) characterized by the prop- 
erty \I/(x; A) = P(x; A) cxpVXx with a polynomial in A function P (see also 
Darboux (1915: p. 212)); 

6. The time isospectral deformations of these reflectionless potentials being gov- 
erned by KdV-dynamics (Ablowitz & Segur 1981); 

7. Lame's potentials u = n (n + 1) p(x) (see Whittaker & Watson (1927)); 

8. Generalizations of the Lame potentials (Darboux 1882, 1915: p. 228) 

, ,.,nCn 2 x ... ,, 9 v{y— 1) ... . dn 2 x 
u = H(H- 1 k 2 —j- + - 1 k sn 2 x + — '- + v\v' - 1 — =-, 
dn x sn^x cn^x 

being called now the Trcibich-Verdier potentials. The following year after 
Darboux, two comprehensive memoires by Sparre (1883) appeared on further 
generalizations. 

A considerable result of Matveev and Its is that they all, apart from the cases 
2-3, lie within the framework of unified theory and are specifications or degenera- 
tions of a wider class, the class of algebraic-geometric potentials expressible by the 
following nice formula (Its & Matveev 1975a, b; Matveev 1976) 

d 2 

u = -2 — \n<d(xU + D) + const. (2.1) 
dx J v ' 

The same authors (see also Novikov 1974; Dubrovin 1975; Lax 1975) proved that all 
these potentials are solutions of higher stationary ode's which are presently named 
Novikov's equations. An important observation of Novikov (1974) was that these 
equations are representable as Hamiltonian finite-dimensional dynamical systems 
in x. Soon this property was completely clarified by Gel'fand & Dikii (1975, 1979). 

The Hamiltonian treatment of Novikov's equations closely joined the examples 
mentioned above with the well-known Liouville's integrability. Moreover, Liouvillc's 
integrability of these equations has received a natural completion. Namely, all the 
solutions are given by the famous trace formula (Matveev 1976) 

g 2 3 +i 

« = 2j> fc -£lSfc- (2.2) 

k=l k=l 

Notice that presence of Liouvillc's attributes itself does not automatically provide 
a procedure of the integration (this is a theorem of existence) but the trace formula 
(2.2) supplemented with the Jacobi inversion problem brings about such a proce- 
dure. In the language of dynamical systems this means, in fact, transformation to 
separability variables. 

What could one say about solution $ corresponding to the potentials mentioned 
above? The answers are well-known. Elementary functions in the cases 1 and 4-6, 
Airy's special functions in the case 3, the Hermite polynomials (up to an exponential 
factor) for the case 2 under A = — 2n — 1 with integral n, Weber's (or parabolic's 
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cylinder) special functions for this case under arbitrary A, and elliptic and related 
to them functions for the cases 7-8. For example, in the case 7 with n — 1 we have 

nx-A) = a{X+ a ^ iX)) e^- 1 ^. (2.3) 

The 0-representation of the ^P-function corresponding to the potentials (2.1), in 
traditional notations and terminology, is given by the formula of A. Its: 

*(^(P))=»„s t (P)5(^±^) e o m , M 

On the other hand, the elementary solutions to the are representable by ele- 
mentary functions, i. e. exactly solvable in terms of elementary indefinite integrals. 
The elliptic cases 7-8 and (2.3) are not exceptions: they are representable by indef- 
inite elliptic integrals. A natural question arises: what about (2.1) and (2.4)? The 
answer is that these are not exception as well. Indefinite quadratures is a common 
property of the spectral problem under consideration. This means that the known 
Liouville's integrability of nonlinear Novikov's equations, in fact, turns out to be 
equivalent to the quadrature integrability of the 'linear' ^ and conversely. In turn, 
the well-known effective solvability of direct /inverse spectral problems in the class 
of analytic nonsingular decaying soliton potentials (Ablowitz & Segur 1981) turns 
out to be nothing else but explicit solvability by elementary tools, i. e. indefinite 
integrations and their inversions. Moreover, such an opportunity is only one. 

The next sections contain proof of the statements above and we suggest that this 
is a common feature of all spectral problems arising in the soliton theory. It should 
be emphasized here that the representation (2.4) (after crossing out A and adding 
the KP- variables y,t 1 ,t 2 , . . .) as a function of a point V on arbitrary algebraic curve 
is a natural and fundamental object for KP-hicrarchics of (2+l)-PDE's (Krichever 
1977). We shall restrict our consideration only to spectral problems defined by ode's 
as independent objects, so that their integrability, in the above sense, belongs to 
their intrinsic nature. Moreover, we will restrict ourselves only to the Schrodingcr 
equation (1.1) that we view merely as a differential equation with a parameter 
(subject to explicit integration) rather than a spectral problem with any boundary 
conditions, commuting operator, etc. 



3. Jules Drach 

Perhaps the most surprising facet is the fact that the ideology mentioned above has 
not received mention in the modern literature in the context. It belongs to J. Drach 
(1919) and his name was revealed by D. & G. Chudnovsky (1984) and V. Matveev 
(Bclokolos et all 1994: pp. 84-85), who drew attention his remarkable results. The 
first sentences in (Drach 1919) clearly indicate his motivationsf . 

f . . . where h is an arbitrary parameter. The most interesting among them are those where the 
Riccati equation 

p' + p 2 = <p + h 

d 2 

(and consequently the equation = [<p{x] + h] y as well) can be integrated by quadratures; we 
will show how one determines the function ip in all these cases. 
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Theorem 3.1 (Drach's Thesis). The class of finite-gap (algebraic-geometric) 
potentials is the only one when the spectral problem (1.1) is integrable for the W by 
quadratures under all values of the parameter A. The solution has the form 

X 

*±(x;A) = y/R^ expj^^, (3.1) 
wherein the function R is a polynomial in A solution of the equation 

A* 2 = -^RR" + ^R' 2 + {u + X)R 2 . (3.2) 

Drach himself does not give a proof therefore to produce one, including step- 
by-step mechanism of integration, is not without interest. The proof will call to 
mind the classical result of Its & Matveev (1975) about criteria of the potential to 
be finite-gap. However we do not invoke any attributes of spectral theories: reality 
of the potential, periodicity, Weyl's bases, monodromies, squares of eigenfunctions 
R = + b^{^> 2 + cV ^2: spectrums, resolvent, functional spaces, etc. 

Actually we will be writing known and somewhat new formulas of the theory but 
keeping in mind only ideology of Drach. Before passing to the proof of the theorem 
we will formulate a known statement which occurs in the literature in numerous 
contexts (Ermakov 1880; Marchenko 1974; Its & Matveev 1975; Gel'fand & Dikii 
1975; Al'ber 1979). 

Proposition 3.2. Integrability of the equation (1.1) is equivalent to integrability of 
the third order linear differential equation 

R'" -4(u + X)R' -2u'R = (3.3) 

or compact nonlinear equation of Ermakov for quasi- amplitude E — >/~R of the <3/: 

H"-(u + A)S = -^. (3.4) 

Though this proposition is a straight consequence of the equations (3.1-2), we 
will give a derivation all the formulas (3.1-4). Arguments are as follows. 

One sufficient test of quadrature integrability a given ODE is furnished by a 
solvable Lie point symmetry of the equation. This technique (group analysis of 
differential equations) is rather well developed (Ibragimov 1985; Eiscnhart 1933) 
so that applying this simple computations to the Schrodinger equation (1.1) we 
get necessary attributes of the theory in a natural way. Indeed, generator <J5 of the 
point symmetries (x, A (x, \&) for the equation (1.1) has the following form: 

8 = (a* + 2 R) d x + {a x V 2 + (R' + const.) * + b} d 9 , (3.5) 

wherein a = a(x;X), b = b(x;X) are arbitrary solutions of (1.1) and R = R(x;X) 
satisfies the equation (3.3). The functions a(x;A), b(x; A) do not help in further 
integrating (need to know solution itself) and we set them equal to zero. The 
remaining free constant in (3.5) says that the symmetry does not disappear and we 
have a solvable commutative 2-parametric symmetry ( <8» 1 — \& , C? 2 — 2Rd x + 
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R'tydy). Transformation to 'integrable' variables (x, ty) (z,w) is provided by 
the standard Lie symmetry machinery. We get an explicit change (Eisenhart 1933: 
p. 91, case 3°) 



X 




In new variables, as Lie's theory guarantees, the equation will be easily solvable. 
Making use of one-fold integrated form of the equation (3.3), i. e. the equations (3.2) 
or (3.4), an equation for the function w = w(z) becomes w zz + w 2 — /i 2 , where the 
constant of Ermakov-Drach fi is a constant of integration. The equation is readily 
solved indeed and, after back transformations, we arrive at the formula (3.1). Point 
of departure of Ermakov (1880) was the same — integrability by quadratures. 



4. Integrability of the Schrodinger equation: Theorem 3.1 

(a) Proof 

Necessity. The function "J" followed by the R is a function of two variables x 
and A. Our aim is to find out the dependence of these functions upon both the 
variables. They are analytical entire functions of the parameter A (Yakubovich & 
Starzhinskii 1975). This allows one to represent R in a form of analytical series 

R{x;X) = Rq(x) + Ri(x) A + R 2 {x) A 2 + ■ ■ ■ + R k {x) X k + ■■■ . (4.1) 

The equation (3.3) is linear hence there exists a recurrence relation on the coeffi- 
cients R k derivable from the equation (3.3) 

X 

1 If 

R k _ 1 =KR k , K = -d xx - u+ - u x . .. dx . (4.2) 

The next step is to make use of Ermakov-Drach's equation (3.2) wherein the 
constant ji is, at the moment, arbitrary and independent of A. We must require 
integrability under all values A. This means that after substitution the series (4.1) 
into (3.2), coefficients in front of A fe must be zeroes, independently of each other: 

( M 2 + \R Q R{ - \R' Q 2 - uRl) + (\(RoRi)" - \KK -2uR R x - i? 2 )A + • • • . 

Hence, the coefficients R k (x) are determined by subsequent integration this infinite 
system of chained equations. It is clear that quadrature integrability takes place 
not for arbitrary functions u — <fi(x). By which restrictions are such functions 
distinguished from all possible ones? 

Rewrite the series (4.1) as a power series in ( = A -1 (Gel'fand & Dikii 1975): 

R(x; A) = Ro(x) + Ri{x) ( + R 2 (x) ( 2 + ■ ■ ■ + R k (x) ( k + ■ ■ ■ , 
whereupon the recurrence relation (4.2) acquires the well-known computable form: 

Rk = & Rk-i ■ R(x; A) = 1 - Q u - c^j C - Q u xx - ^ u 2 + i c x u - c 2 ) c 2 + • • • • 
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Thus, for arbitrary u = 4>{x) the \P is expressed by the formula (3.1) with the differ- 
ential polynomial R([u]; A) of infinite order. Moreover, making use of this recurrence 
with subsequent collection in (, the equation (3.2) turns into 

C^ 2 = 1 + 2 Cl C+ (2c 2 + cf)C 2 + 2 (c 3 + c lC2 )C 3 H h second half. (4.3) 

The dependence on u has gone at infinity: the first half (it is infinite!) contains 
only constants c k and the second one does differential polynomials of infinite order 
in u. In order to get conditions of finite order in derivatives we have the only 
possibility — to set R to be a finite polynomial: 

R([u]; X) = X 9 + Rg-^- 1 + --- + R 1 X + R a . (4.4) 

Accordingly, splitted restrictions on the potential (4.3) become g, now not necessary 
zero, constants I k = F k (c l7 c 2 , . . . , c g ) plus g+1 differential conditions of finite order 

I k =F k (u, u x ,..., u (2s) ; c l7 c 2 ,..., Cg ), k = g + l...2g+l. (4.5) 

Their compatibility means that such a class of potentials is not empty and the 
parameter fi must depend on the spectral one. This corresponds to a choice the 
particular solution R among three linear independent -Ri,2,3 defined by integration 
constants ^1,2,3 such that the constant ^(^2,3) is algebraically related with A: 
/i = /i(A). These potentials admits such a dependence and the condition (4.3), or 
which is the same (3.2), turns into a hyperelliptic algebraic curve of finite genus 

<&: ^ = X^+hX^ + --- + I 2g X + I 2g+l = 
= (X-E 1 )-..(X-E 2g+1 ). 

The set of equations (4.5) is 1-fold integrated Novikov's equations (stationary 
KdV-equations) and a half of their integrals I k = F k ([u];c). Thus quadrature in- 
tegrability of the ^ has became equivalent to integrability of these equations. If it 
has been established, a final answer is given by substitution u = 4>{x) into the poly- 
nomial (4.4) and then into the formula (3.1). We emphasize that the well-known 
formula (3.1) does not mean any integrability without formula for (f)(x). This is 
ansatz and such a form of solution to the ODE (1.1) does exist for arbitrary <f>(x). 

Sufficiency. The sufficiency is a procedure of explicit integration of the equations 
(4.5). It is widely known and we do not repeat details here (Drach 1919; Dubrovin 
1975). The answer is the formula (2.2) plus integral representation for the functions 
7 k {x). Nevertheless we note that being a remarkable identity in the spectral theories 
(Levitan 1987), the trace formula (2.2) turns into a necessary and key object in the 
quadrature methodologyf . □ 

(b) Drach-Dubrovin equations and formula for the ^-function 

Rather than mere identities, we consider consequences of the well-known definition 
of the polynomial R and new variables 7 fe , for example u xx — 3 u 2 = 7j7fc> as al- 
gebraic transformations from the variables {u, u x , . . .} to the separability variables 

f Independent role of these formulas for the KdV-cquation, as well as their dynamical relatives, 
was pointed out by Matveev (1975). See also his appendix to (Dubrovin et all 1976). 
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{7, /i}. Namely, writing the function 
i?( M ;A)=A«-(|- Cl )A«- 1 -(i 



(4.7) 



in factorized form R([u]; A) = (A — -f 1 (x)) ■ ■ ■ (A — 7 9 (x)) , we rfo define the first half 
of the change of variables as zeroes of this R: 



{u,u X7 ...,u {2g 1) } i-> {( 7l , 7 2 , . . . ,7 fl ), (^,^2,...,^)}. 
More precisely, the first part of the change (4.8) is as follows 



(4.8) 



29+1 



fc=i 

3 



fe=l 



u ra = 12 ^ 7 j? - 8 ^ 7fe7j ~ 16 Cl ^ 7fe + 4c? + 8c 2 



fc=i 



fc=l 



and not closed due to the odd derivatives of u. Missing ones and therefore the 
second half is extracted by involving the Drach-Dubrovin differential equations: 



-2 H 



d 7 fc = 

Ax II(Tfc-Tj)' 



Ai?) - (7 fc («) - ^1) ■ ■ ■ (7fcW - ^2 fl+ i), (4-9) 



where constants Ej are functions of Novikov's constants Cj and the integrals (4.5). 
We thus get the remaining part of the complete change (4.8): 



-4E 



fr[ II (Tfc-Ti) 



i 2£7 fe + 2c 1 



fe=l 



^ n (Tfc-7,0 



fcJ>fc 1 n (t*-t») n -t«) 



The equations (4.9) are readily rewritten into a promised integral form which co- 
incides the Jacobi inversion problem for the hypcrclliptic algebraic curve (4.6) 



(7i>Mi) 



/ 



dz 



(7g, Mg) 



+ ... + 



/ 



(7i. Mi) 



dz 
z — + 

w 



(7g,Mg) 



+ 



dz 

10 

dz 



= di 



(4.10) 



(7i> Mi) 



(7g> Mg) 



_! dz 



dg — 2x 



Article submitted to Royal Society 



What integrability means? 



9 



With this equation we can obtain formula for the ^-function (3.1) which is not 
seem to be in the literature: 

The variables of integration lie on the curve w 2 — (z — E\) ■ ■ ■ (z — £?2 S +i). 

It may be remarked at once that complete set of the transformations above to the 
separability variables {7, fi} is not necessary for Novikov's equations themselves. 
This is an attribute of their Hamiltonian description (Gel'fand & Dikii 1979) deriv- 
able from the main trace formula (2.2) and the equations (4.9). It thus appears that 
the formula (3.1) or (4.11) supplemented with the two objects (2.2) and (4.9-10) 
explains the nature of integrability in question. Subsequent procedure of inversion 
for a symmetrical sum of 7's is necessary for ultimate representation to the solution. 
In the degenerated cases (constant potential, solitons and the like) all the integrals 
(4.10) reduce to integrals of rational functions, so that the problem becomes trivial 
(inversion of logarithms) and leads to well-known exponents. In general case, this 
transcendental problem is solved with the help of Riemann's 6-functions (see § 7) . 

We would like to note that even 'naive' (not finite-gap) procedure of integration 
of the trivial potential ^" — const. ^ = X^ completely accords with the scheme 
described above. Holomorphic integrals do not appear but, instead, the following 
objects arise: integral definition of the logarithm, i. c. integration of a trivial rational 
function, and necessity of inversion of the former (the exponent). 

5. Some consequences 

(a) Integrable X-pencils 

Quadrature integrability docs not depend on a choice of dependent/independent 
variables. We could formally avoid the inversion of Abelian integrals (4.10) rewriting 
the theory in the 'inverse' variable u. Indeed, after the change (x, "J) <-> (u, Y) 

x = X {u), ^ = ^\—Y (5.1) 

we arrive at a second order linear ode ({%, u} is the standard schwarzian) 

Y uu = Q(u;X)Y, Q(u-X) = - l -{x,u} + {u + X) x 2 u (5.2) 

which can be of interest in its own right as a new spectral problem (operator A- 
pencil) for some 'good' functions Q(u; A). The former depends on a chosen potential 
x = x( u )- Apparently the second equation in (5.2), written in the form 

FF UU -^F 2 -2(u + X)F 4 + 2Q(u)F 2 = 0, F = Xu (5.3) 

and equation (6.4) play an important part in the theoryf , if only because inversion 
of arbitrary finite-gap potential satisfies this integrable 2-nd order nonlinear ode 
with suitable function Q. Let us consider some examples (JV = n (n + 1)). 

f See also a footnote on p. 97 in (Gel'fand & Dikii 1975) apart from two misprints in one term. 
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• Soliton potentials u = —n (n + 1) cosh 2 x: 



3 1 A-l X-N-l 

Q(u;X) = - — 



I6(u + N) 2 Au 2 Au{u + Ny 
• The Lame potentials u = n (n + 1) p(x; g 2 , g 3 ) (Whittaker & Watson 1927): 

1 (3 + 2N)u + 2NX 



16 ^ {u _ ^2 8 (u _ Nei){u _ iv e2 )( u - 7Ve 3 ) ' 

• Arbitrary even elliptic finite-gap potential u = U(p(x)) , where U is a rational 
function of p. A wide family of such potentials there provides the theory of 
elliptic solitons (Acta Appl. Math. 1994). The function Q is as follows 

_ 1 {U, p} ^ U(p(x))-3p(2x) + X ~ 
Q=^ ± J7T ± + 77T- 2 =U(p{x);X), 

2 u p U p p 2 x 

where U is another rational function of p. Hence Q = J7(J7~ 1 (u); A) is a genus 
zero algebraic function of u or rational function of p: Q = Q(p; A); 

• Arbitrary elliptic soliton u = <fi(x). The equation (5.2) takes the form 

<&uv ®u 1 + $vv $u 1 U + X' 



V<$>1 2 w $3 4 W 2$2 



where ^(u, v) = denotes a differential equation connecting the elliptic func- 
tion it and its derivative u x = v (algebraic equation of genus unity). 

By the previous constructions, all these equations and their relatives of the type 

(3.3) are of Fuchsian class and integrable by quadratures. We thus get a solution for 
all elliptic solitons generalizing formulas of Hermite (Belokolos et all 1994: p. 82): 

u 

y( U ;A).^K^exp/^^, (5.4) 

wherein, owing to (4.7), 1Z(u, v;X) = v R([u]; A) becomes rational function in (u, v). 
The last step is a standard problem to representation the elliptic Abelian integral 

(5.4) , which is solved in terms of Jacobian ^-functions. 

It should be remarked that the presence of explicit <J/ leads to explicit factoriza- 
tion all of the linear operators (1.1), (3.3) and (5.2) with arbitrary A's. For example: 

d xx - (u + A) = (d x +p)(d x -p) = 0, 

where roots ±p(x; A) have the quadrature form 



1 B' 4- „ 9 1 9 



The 0-functional representation for this and other factorizations is readily written 
down using formulas of § 7. 
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(b) Liouvillian integrabilities 

To all appearances Liouville was the first (1833-41) to recognize a significance of 
the object R = ^f^ 2 an d associated linear ode of higher order in the context of 
integration of linear ode's in closed form (Liouville 1839). The third order equation 
(3.3) explicitly arose and was discussed on pp. 430-431 in (Liouville 1839). Though 
Liouville was doing in the spirit of algebraic solvability!, a presence of a parameter 
in equations turns the theory into the spectral one. The polynomials in A arose 
also in works by Darboux and polynomial in p(x) was considered by Hermite in 
the case of Lame's potentials u — n(n + 1) p(x) (Whittaker & Watson 1927). 
This corresponds exactly to the 'polynomial in it' cases of Liouville (1839) for the 
equation (3.3). 

V. Kuznetsov (2001, personal communication) pointed out a relationship of the 
theory with an algorithm of Kovacic (1986) and, as we have seen now, this key 
observation leads to the natural conclusion: 

• Liouvillian integrability of linear ode's (1830-40's) with a parameter in fi- 
nite terms is equivalent to algebraic integrability by Liouville (1840-50's) of 
nonlinear Hamiltonian systems]. 

The theory and examples above show that this is not a coincidence and, probably, 
the modern efficient computer-algorithmic theory, being applied to both equations 
(1.1) and (3.3), would provide the independent approaches to generation/classification 
of intcgrable linear operator pencils. For excellent explanation of Liouville's ideas 
see books by Mordukhai-Boltovskoi (1910), Ritt (1948), references in works by the 
authors mentioned in the last footnote and the Chapter IX in (Liitzen 1990). Among 
other things this book provides a full account of references for further study. The 
next section contains an additional information, examples, and connections with 
non-intcgrable equations. 

6. Integrable cases of Riccati and related equations 

For a simplicity and to avoid lengthening the terminology we will refer to the second 
order linear ode's (or potential) and corresponding to them Riccati's equations of 
general form y z +a(z) y + b(z) y 2 = c(z) as one object. Well-known transformations 
between them have a quadrature characterization. 

(a) Riccati's equations 

The Q-functions corresponding to integrable Riccati's equations (5.2) can be ratio- 
nal/algebraic, elementary, or transcendental. It is rather evident that non-integrable 
equations contain integrable subcases. Say, the quantum harmonic oscillator Y uu = 

f See however p. 456 in (Liouville 1839) about what is nowadays named Liouville's extension. 

f We speak here only about general link and avoid discussion the rigorous correspondence 
between finite-gap operators and Liouvillian solutions, differential algebra (Ritt 1948), algorithms 
of Singer (1981), Kovacic (1986), Picard-Vessiot theory (van der Put & Singer 2003), etc. In 
particular, we do not touch an important question: when does isomorphism between these two 
Liouvillian integrabilities take place? We should mention here some comments about this analogy 
in (Morales-Ruiz 1999: pp. 51-52) (see also van der Put & Singer 2003). However main attributes 
of the theory (spectral curves, polynomial in A, ©-functions, etc.) are not discussed in these works. 
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(?i 2 — 2/i — l)Y shows that this potential, seemingly having nothing in common 
with the finite-gap ones, is intcgrable by quadratures if /x is an integer. 

Let us reverse a view on the equations and change (5.1-2). Whether exists 
transformation between nonfinite-gap equation (with a parameter or no) 

Y„=Q(z)Y (6.1) 

and a finite-gap one? Such a transformation depends on chosen equations and can be 
rather complicated. In contrast to the preceding (5.1-2), corresponding functional 
relation is given by the general change of variables (x, \I>) (z, Y) 

x = X (z), * = Vx~z Y (6-2) 
and depends on the potential u = (f>(x)- 

Proposition 6.1. Arbitrary equation (6.1) and the finite-gap one (1.1) are trans- 
formable into each other by the following functional relation 

X - * 2 (*;A) Y 2 (zY [b - 6) 
where functions ^ 12 {x;X), Y 12 (z) are independent solutions of (1.1) and (6.1). 
Proof. From (1.1) and (6.1-2) we have 

~{X, *} + (*(*) + A) x2 = <?(*)■ ( 6 - 4 ) 

Clearly, the sought for functional relation is an integral of this equation. The po- 
tential <j>(x) is defined by the corresponding ^ which is known. From the second 
equality in (6.2) we have 

dx = dz 

Integrating and supplementing with the property 

dx * 2 

we get the formula (6.3) and complete the proof. □ 

It thus appears that the product solution R = * 1 ^ , 2 1S a fundamental object 
in the finite-gap theory and the ratio (6.3) is fundamental in transformations be- 
tween Riccati's equations. Such arguments might seem to be trivial because every 
intcgrable equation is transformable to the trivial Y zz = 0. But an example of Raw- 
son (1883) shows nontrivial consequences: generating of finite-gap spectral problem 
(1.1) with <f>(x) — n (n + 1) x~ 2 from the oldest and classical equation of Riccati- 
Bernoulli with a parameter. The paper of Rawson is so short that we completely 
reproduced it in the Appendix A without any comments!. 

The following extra examples exhibit a functional relation between the Lame 
potentials and equidistance spectrum of harmonic oscillator Y zz = (z 2 — 2 \i — l)Y. 

f See also Liouville's (1841: pp. 11-13) classical considerations on intcgrability of the potential 
<j>(x) = Bx~ 2 with appearance B = n(n + 1). Another example of 'triviality' is the class of 
potentials in elementary functions (solitons and the like) generated from the zero potential by the 
Darboux transformation (Matvccv & Salle 1991). 
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• <f)(x) = n (n + 1) x 2 and \x = 0. We obtain 



z 

X : x 2n+1 = J e z2 dz, A = 0; 



The Lame 1-gap potential <j>(x) — 2 p(x). The relation (6.3) has the form 
cr(a + a;) J H 2 (z) 



As we have seen, these formulas are representable by quadratures in both vari- 
ables. List of examples along these lines may be readily extended. For example one 
to derive a generalization of Rawson's transformation (see Appendix A) from the 
classical equation of Riccati to the Lame one. A counterexample of Liouville-Airy 
Y zz = (z + a)Y brings out the differential field independently of the parameter a: 



f 2^dx f 



dz 



(z + aY 

as might be expected for the variable z. 

In a broad sense, the variables u = (j){x) and &(x; A) may be thought of as 'con- 
venient' variables for all integrable Riccati's equations with a nontrivial parameter 
because, in this case, the differential polynomial R([u]; A) has a universal description 
as u = <t>{x) is a solution of the equations of Novikov. 



(6) Related equations 

The equations of Riccati, Ermakov-Drach, Novikov's equations, and equations of 
the type (5.1-3) and (6.4) are hidden forms of integrability of one another but are 
not only integrable linear/nonlinear equations arising in this theory. Besides them, 
the ^-function itself satisfies some nonlinear homogeneous autonomous differential 
equation of the third order with a parameter(s). That equation and its relatives 
are obtained with the help of suitable elimination of the potential u. The following 
instances illustrate the remark above. 

• <p(x) = n (n + 1) p(x + c): 

n (n + 1)(* - *'*") 2 = 4 *(*" - A x *) (*" - A 2 *) - A 3 *), 
where arbitrary parameters \ k are restricted by the relation X 1 + A 2 + A 3 = 0. 

• Arbitrary elliptic finite-gap potential u. Then the \1/ satisfies the equation 

where $(u, v) — is an algebraic relation between u and its derivative u x = v. 
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The result of elimination depends on a chosen independent variable x, u, . . . and 
all these integrable equations and their i-deformations can be of interest in their 
own right if only because they are closely related to known 3-rd order autonomous 
nonlinear ode of Jacobi for the ^-constants in the framework of Fuchsian equations. 
We do not develop this topic here. It would appear reasonable that the widely known 
and universal 9-description of the theory (Krichever 1977a, b) should be obtainable 
from the spectral problem itself. This is so indeed. 



Here we will obtain the 6-representation (2.4) for the ^-function and, thereby, its 
properties as a function of Baker-Akhiezer. Insomuch as regular derivation of this 
axiomatic representation is not described in the literature, one to write up that 
procedure is, perhaps, not without interest. 

As we mentioned above, pure spectral approaches were extended to the complex 
valued potentials. (See recent monograph by Gesztesy & Holden (2003) for most 
exhaustive bibliography and new results in spectral treatment of the theory). Taking 
this into account we will refer a general A-dependence of the as its spectral 
property and its dependence upon x-variable is considered to be parametric. Such 
a spectral view arose in the paper by Akhiezer (1961) and was completed in full by 
Its & Matveev (1975: §4). It seems helpful to compare this classical approach with 
a 'reverse' one, i. e. primary ^-dependence of the \&. This would correspond to pure 
quadrature arguments with a parametrical A-dependence made out in the previous 
sections. 

Theorem 7.1. The 6 -functional representation (2.4) to the \P -function is a con- 
sequence of the quadrature representations (3.1) and (4.11). 

Proof. Since the spectral parameter A is connected with the variable /j, by the alge- 
braic equation (4.6) of finite genus g we view both these variables as meromorphic 
functions A = A(r), fj, = fi{r) of a global parameter r on the curve (4.6). Accord- 
ingly, we consider the ^-function (4.11) as a function of x and r 



where a k arc arbitrary constants. It is a symmetrical function of the quantities 7^ 
and the formers, as functions of x, are defined from the inversion problem (4.10): 



7. 




7fc0*0 



(7.1) 




1 



(7.2) 



a. 
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Let Ajk be a matrix of a-periods of the holomorphic Abelian integrals (7.2): 

)du j (z) = A jk . 



All the terminology and notation in this section is standard and elucidated in any 
paper on the finite-gap integration. As usual, we normalize the integrals Uj(z), 
introducing the canonical base of normal holomorphic Abelian integrals oi{z): 

Wj(z) = A jk Z k (z), (j)dti j (z) = 8 jk , j)&u {z) = B 3k . (7.3) 

a k b k 

Jacobi's inversion problem (4.10), (7.2) acquires the form 

= -xUj + Cj, Uj = 2(A~ 1 ) jg , (7.4) 

k=i 

where we adopt the concise Baker's notation (Baker 1897) for Abelian integrals: 

tk 



The vector U depends on the curve & but a k and Cj are arbitrary constants. The 
function (7.1) depends on the point r on the curve & and on the variable x through 
7's. Thus it is a single- valued function (but not Abelian) of a point r)(x) = —xU+C 
on the Jacobian Jac(&) and, hence, has a 0-functional representation. Indeed, the 
arisen sum (7.1) is nothing else but the fundamental T-function of Weierstrass 
(1856) and Clebsch-Gordan (Clebsch & Gordan 1866) 

7i 7 2 7 S 

tJ^ ■ ■ ■ > 7 9 \ = r dfiej)(z) + r d ^ (z) + ...+ /" d n ?7)( ^ 5 

\a l7 a 2 , • • • , otgj J J J 



where dll^^ (z) denotes the elementary normal Abelian differential of 3-rd kind with 
1-st order poles at the points z — £, f] and residues +1,-1 respectively. Clebsch and 
Gordan (1866) devoted Chs. 6-8 in their book to detail properties of this object and 
regular procedure of derivation the 0-representation for it (see also Baker 1897). 
We have from there (Clebsch & Gordan 1866: §§54, 57; Baker 1897: §§171, 187-188) 



IL x a 'l = U a ^{x)-U a ^{z) = ln 



9(^'"+r) / e(S''°+r) 



yi,m 1 _ . . . _ ~zg-i,m s -i j g a zerQ Q £ ^ e 0-function. 
In our situation we have a = A, v = oo, x = j k , z = a k . The representations for 
zeroes r and the Riemann constants 7C are not unique because they depend on 
a lower bound of the holomorphic integrals u>(z). On the other hand the point r] 
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on the Jacobian contains free constants Cj in (7.4), so that we may simplify the 
considerations putting m's equal to 7-'s apart from j k and set 

r j =ui j {V 1 ) + ---+u; j {P g _ 1 )+lC j 

with arbitrary points Vj, say 7's. We thus arrived at Riemann's function <d(u>(P) — 
e). As a consequence of these arguments and explicit zeroes of the >J> in (3.1) wc 
get the following formula (see also Baker (1928): pp. 588-589) 

n A. - + ■ ■ ■ + n A : = m e(s(7i) + ... + g(7g) _ g(00) + x:) + 

The integrals in (7.1) are elementary but not normal. Hence we have an identity 

g 7fc ^ g 

with some normalizing constants h- k (r). Taking into account Jacobi's problem (7.4) 
we arrive at the intermediate answer 

, N e(u(X(T)) + xU-C-K) » , / \ \ 

- e(a(c,) + ^-c-x) <^C V'HUM) ■ 

Meromorphic part of the ^-function has been determined. In order to determine 
the function f(x;r) = ^""^ (r) u)j (j k (a;)) in (7.5) wc involve the transformation 

properties of the IP as the function (7.1) on the curve &. Let a k r and 6 fe r de- 
note linear-fractional transformations of the variable r corresponding to the cycles 
a k and b k . From the primary representation (3.1) we conclude that ^/(x;t) saves 
its own form (3.1), up to a multiplier M(t), when r undergoes the transforma- 
tions 2, b. Hence the formula (7.5) must hold this property. Clearly, the function 
f(x;r) is an entire function of x since uj(j)'s are everywhere finite. Further, the 
^-transformations of the O's in (7.5) say that the function / has to be a linear 
function in x to compensate an exponential multiplier c(r) exp(— 2-7ria; Uk) in (7.5): 

f(x; t) = k(t) x + const(r) . 

Invariance of the 9-functions in (7.5) with respect to 2-transformations implies 

1): x{a k T) = x{T). 

The transformations b imply 

2): x(b k r) =x(T)+2niU k . 

Let t approaches the pole of the meromorphic function A(r). Such a pole is only 
one and, designating £ = r — , we have 

A 2 B . . (A\ 2g+1 2 3 + l B ( A^ 2g 
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From these formulas and (3.1) we obtain (extracting terms independent of 7's) 

±/i da; 



ln*(x;T) = i^ln|^ + | + (C- 7fc ) + -"} + 



(A- 7l )...(A- 7s ) 

— W-{^ + -} ± /{S?« + -}d. ± {f + ^ + ...}. 

and therefore we get one more property 

3): ±«(r). 7 A_ + »+.... (7.6) 

A function with the properties 1-3) does exist on &. This is the normal ele- 
mentary Abelian integral of the second kind ±x(r) = O(r) with the first order 
pole at the point t x . Its principal part and a constant in (7.6) are well defined so 
that this function is completely determined and unique. In turn, the three param- 
eters {t^ A, B} may be freely chosen, say, {0, 1, 0} respectively. The vector U 
in (7.4-5) becomes the vector of b-periods of the integral f2(r). Summarizing the 
arguments above and recovering a normalizing constant 

*(*■ a(t)) = e (*M- g ) e(*{KT)) + *u-D) e n {T)x (7 7) 

[ ' 1 )] e(«(A(r))-D) 6 (£(00) 1 j 
we get the spectral properties of the \& and complete the proof. □ 
From the last expansion of the ln^x; r) we obtain the formula (2.1). Indeed, 

d /u., T ,^ o^.\l _ d 1 e(«(A(r)) +**/-£>) 



ff , dr . . o/ \ II d , 6wAr +z£7- 



hence, by virtue of the property [/ — 



dr 

replaced by —djdx and we get the final formula 



the derivative 4r 



may be 



, 2 2g+l 

u = -2— ^ lne(a;[/ + £(oo) -23) - ^ £ fe . (7.8) 
dx fc=i 

We should conclude here that the spectral and quadrature considerations are 
mutually replaceable. The explicit transition between these approaches there pro- 
vides the known Weierstrass's theorem on a permutation of arguments and param- 
eters in the normal Abelian integrals of 3-rd kind 

fL^Oz) - ^v.A x ) = n«, x (^) - n z>x (/i). 

This fact immediately leads to the equivalence of the 'spectral' formula for the ^ 
(Its & Matveev 1975: formulas (4.12-13)) and the quadrature one (4.11) or (7.1). 
The above mentioned theorem, in our notation, has the form 

1 f w + (i dz 1 f ( w + ji k (x) w + [3 k \ dz 



f ( w + fi k {x) w + f3 k } dz 

/ S 7 's - 7 > h holomorphic part (a;), 

J K*-7 fc (a0) ( z - a JJ w 



2 7 (z-A) u; 2 7 l(z- 7fc (aO) 
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where (5\ — (a k — E\) ■ ■ ■ (a k — E2 g +i). Since the contains all the information, 
we can obtain suitable expressions for all objects of the theory. In particular, 0- 
functional representation for 'finite-gap function' of Ermakov (3.4) having numerous 
applications. Renormalizing it by the formula S 2 = ^+(x; A) A) we obtain 



'Finite-gap' means that the constant \i in Ermakov's equation (3.4) is not indepen- 
dent of A: (i = /i(A). Note that the equation (3.3) itself exemplifies the integrable 
and factorizable linear operator pencil with Abelian coefficients (like soliton spectral 
problems) but its solution, as a counterexample, is not a function of Baker- Akhiezer. 

Many of constructions in §§ 3, 5-7 can be carried over to arbitrary spectral 
problems although not so simply as in the case of the Schrodingcr/KdV equation. 
Nontrivial examples of the \I/ were obtained in (Ustinov 2002) and modifications of 
Dubrovin's equations and trace formulas in (Brezhnev 2002). 

8. Conclusive comments and bibliographical remarks 

Apparently Baker (1928: p. 587) was the first to realize the exponential property of 
the fundamental T-function in disguise, however explicit meromorphic integrals like 
f2(r) are, to all appearances, the result of the modern theory and became a universal 
property of all the integrable models (Krichever 1970's). Note, that throughout the 
paper we made no restrictions on the curve (4.6). It may be singular and the theory 
can be rewritten (with minor changes in § 7) because the degenerated holomorphic 
integrals u>, in this case, turn into the integrals II of the 3-rd kind. Clebsch & 
Gordan (1866: §43) call the corresponding problem 'extended inversion problem' 
which is solved by theta as well. Such cases of degenerations were considered also 
by Baker (1897) and even by Abel (see his CEuvres 1881, I: p. 170-. . . ). 



It is to be noticed, that contrary to commonly accepted (?) viewpoint, the general 
©-series (B 1), as well as its argument z = u>(V) — e, ought to be considered not 
as a special function or formal generalization of Jacobi's ^-function but regularly 
derivable fundamental object for explicit representation of all Abelian integrals, 
meromorphic/uniformizing functions, exponential Baker-Akhiezer-functions, and 
the theory as a wholef . Main ideology belongs to K. Weierstrass (1902: pp. 513-538) 
and was expounded in the books by Tikhomandritskii (1885; 1895: pp. 199-232) 
before publication of Weierstrassian lectures on Abelian transcendents, but with 
explicit use of the Riemann surface. Primarily and naturally arising non-canonical 
form of the ©-series is derived in these works with basic properties and identities. 
The canonical formula (B 1), is obtained after the normalization (7.3). See also the 
book by Clebsch & Gordan (1866: pp. 193-198). For lack of space, we do not pursue 
these important points here. This will be written up elsewhere. By the same reason 
we have somewhat reduced exposition between formulas (7.4) and (7.5). 

f Jacobi's case g = 1 is not exception (Tichomandritzky, M. 1884 Math. Ann. XXV, 197-202). 
Riemann and Baker (1897) do not elucidate an origin of appearance the 0-series. Hypcrclliptic 
case, based on Weierstrassian lectures, was considered in dissertation by Tikhomandritskii (1885). 



S 2 (x; A) = const. 




(a) On the -series 
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(6) Algebraically integrable Hamiltonian systems 



In what concerns an inverse transition in the thesis on p. 11, i. e. transition from an 
algebraically integrable Hamiltonian system to some spectral problem, the answer 
is positive though somewhat noneffective. We assume that separability variables 
{7, fi} sitting on a curve, say (4.9), exist. We can construct the Baker-Akhiezer 
function (7.7) and the Schrodinger operator (it certainly exists) for it. The compu- 
tational part of the Liouville theorem is based on the Hamilton-Jacobi theory of 
canonical transformations. The variables 7's are always taken as poles of the ^: 



since they solve the problem (4.10), (7.4). Ostrogradskii's variables {u, u x , . . .} are 
constructed explicitly by the trace formulas. Isospectral i-deformations of finite-gap 
potentials are readily included in (7.2): the quantities dk become linear functions of 
xU + tV. Accordingly, there is not essential difference between x- and i-equations 
for the ^. Both of these equations (Lax' pair) may be thought of as quadrature 
integrable (factorizable) spectral problems/pencils: Ermakov-Drach's constant /1 
becomes an eigenvalue of the second commuting operator connected through the 
curve with the first one: 



There are infinitely many Abelian functions built from {7, /i} (followed by physical 
coordinates {^(7, fJ,), q("f, /u)}) and expressible through the 6's. Non-constructiveness 
can appear because transformations between various dynamical systems are not nec- 
essary to be canonical. However canonicity is not a necessary attribute of quadrature 
integrability. One fundamental Abelian function determines all the transformations — 
the potential u. 

The author thanks Professor J. C. Eilbeck and Professor E. Previato for numerous discus- 
sions and EP for hospitality in Boston University where the work was carried out. 
The project was supported by NSF/NATO-grant DGE-0209549. 

Appendix A. Note on a transformation of Riccati's equation 

R. Rawson. The Messenger of Mathematics XII (1883), 34-36. 
Riccati's equation 



Jki x ) : 



0(w(oo) + xU - D) =0 



A([u];d x )^(x;X)= f i(X)^(x;X). 




(1) 



is readily transformed into 





by means of the two equations 



= P, a function of x 



(3), 
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If, therefore, either (1) or (2) can be integrated, then the other can be integrated 
also by means of (3) and (4). 

An interesting case of the above transformation is when 

<j, {xi) = hx n + JL (5)) 

x l 

P = x m (6). 

We then have 

dy + 2 = a_Ml_ a (n +2)m -2 + (4 gig + 1) f 1 (g) 

ax a Aax A 

X — X ^ 

m-i m — 1 

2a 1 mx Vi=ay+— . 

2 x 

In equation (7) let c = 0, and n = — - — — , where p is an integer, then (7) and 
(8) become 

dy x -- a- 

dx 1 



T^ "://: J >>■<■. ■ (9), 



dy 9 a,bm 2 ±2™ ? m 2 — 1 ... 



and therefore since (9) is soluble, so also is (10). 
Equation (10) is made linear by putting 

1 dz 

y = — T> 

a zax 

and we thus find 

d 2 z { , 9 a- n m 2 - 1 ] 

or 

d 2 z 1 ,9 ±2(29+1) 9 9(9 + 1) 1 , N 

^2 =|« 1 6(2g+l) 2 a; ^^- 2 + ^^ \z (11), 

where 

m = 2q+l. 
The equations of transformation then become 

_ 2q+l 
— X , 

2 ai (2g+ l)^ 9 ^ = ay + -. 

x 

The differential equation (11) is of some interest as it includes as a particular 
case the well-known differential equation 



C?X 2 I X 2 



Robert Rawson. 
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Appendix B. Notation to § 7 

The 0-function corresponds to the symmetrical matrix Bjk (7.3) 

Q(z\B) = Q(z 1 ,...,z g \B)= e^< BN ' N > +2vi ^ N ' z \ (Bl) 

where (BN, N) = ^ B jk NjN k and (AT, z) = N x z x H h iV g z ff . Canonical base of 

cycles (dj, bj) on © is chosen according to the intersection scheme: dj o b k = Sjk- 
The normal Abelian integrals of the 2, 3-rd kind and their periods have the form: 

Q = — !- + Ci (t - r^) H , (j)dtt = 0, (j)dtt = 2-KWk, 
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